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Abstract: In the current work we study various models of holographic superconductors 
at low temperature. Generically the zero temperature limit of those models are solitonic 
solution with a zero sized horizon. Here we generalized simple version of those zero tem- 
perature solutions to small but non-zero temperature T. We confine ourselves to cases 
where near horizon geometry is AdS'^. At a non-zero temperature a small horizon would 
form deep inside this AdS'^ which does not disturb the UV physics. The resulting geometry 
may be matched with the zero temperature solution at an intermediate length scale. We 
understand this matching from separation of scales by setting up a perturbative expansion 
in gauge potential. We have a better analytic control in abelian case and quantities may 
be expressed in terms of hypergeometric function. From this we calculate low tempera- 
ture behavior of various quatities like entropy, charge density and specific heat etc. We 
also calculate various energy gaps associated with p-wave holographic superconductor to 
understand the underlying pairing mechanism. The result deviates significantly from the 
corresponding weak coupling BCS counterpart. 
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1 Introduction 

In recent years string theory or, more specifically, gauge-gravity duality has seen interesting 
applications in the field of condensed matter physics. One of the earliest such applications 
is the discovery of a superconductor-like phase transition in AdS with a Reissner-Nordstrom 
black hole and a charged scalar field minimally coupled to a U{1) gauge field [1-3]. Other 
system with similar properties include so-called "p-wave" holographic superconductors with 
non- abelian gauge fields instead of a scalar coupled to U{1) gauge fields [4-6]. In various 
works [7-15], authors have studied related systems in the zero temperature limit. 

Generically the zero temperature solutions turn out to be a solitonic solution with 
a zero sized horizon. The authors find that the effective potential for small gauge field 
fluctuations vanishes near the black hole horizon for the abelian cases. This implies that 
the normal component of the A.C. conductivity never vanishes, even at zero temperature, 
which in turn indicates that the superconductor is gapless [7]. However the corresponding 
effective potential does not vanish at the horizon for the non-abelian case. It is concluded 
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that the holographic non-abehan superconductor does have a finite gap for the relevant 
gauge field fiuctuations [14] . The non-abelian system is an anisotropic system which shows 
different conductivity in different directions. 

The near horizon geometry of those zero temperature solutions are interesting and 
ranges from simple AdS"^ to various complicated and Lifshitz like geometries. Near the 
horizon these geometries could be constructed by an analytic perturbation theory [7, 14]. 
Near horizon values of the scalar field or the appropriate component of the gauge field 
enters as a undetermined parameter in those perturbative expansion. Those parameters 
are determined by a numerical integration to infinity and consequent application of proper 
boundary conditions. 

In this work we generalize simple version of those zero temperature solutions to small 
but non-zero temperature T. It should be noted that a non-zero temperature solution was 
already obtained numerically [3, 7, 15]. However from a purely numerical solution it is 
difficult to conclude about the low temperature analytic behaviour of various quantities. 
Whereas we will be able to calculate the nature of various interesting physical quantities 
analytically. We confine ourselves to cases where near horizon geometry is AdS"^. We 
expect that at non-zero temperature a small horizon would form deep inside this AdS"^. 
We intuitively understand this by separation of scales. As the black hole is situated deep 
inside the AdS^ it does not affect the UV physics. Hence we expect that at an intermediate 
scale (rj) the non-zero temperature solution approaches the zero temp one. Importantly for 
a very small horizon size (ro) the intermediate scale may itself be chosen very small so that 
the zero temperature perturbative method would be valid in the intermediate scale. Here 
we have tq ^ ^ 1. We show that we can set up a perturbative expansion in terms of 
the gauge field Aq which interpolates between the black hole horizon and the intermediate 
scale. From this matching in the intermediate scale we argue that a slight variation of the 
zero temperature numerics may be applied to the non-zero temperature case. 

From our solution we may calculate how entropy, specific heat etc. vanishes near 
zero temperature. We also calculate the various energy gap associated with the systems. 
Especially in the non-abelian case we calculate the various energy gap in the system and 
from their ratio we find some hint of underlying "pairing mechanism" . The ratio deviates 
around 33% from its weak coupling BCS counterpart. 

Our results may be generalized to various cases where near horizon geometry at zero 
temperature deviates from AdS"^ [7-12]. The application part may include calculation 
of various fermionic propagator, calculation of second sound and more interestingly low 
temperature behaviour of non- universality of viscosity entropy ratio etc [16-23]. 

Plan of this paper is as follows. In section 2 we discuss the abelian or s-wave case. In 
section 3 we discuss the non-abelian or p-wave case. 
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2 Abelian holographic superconductors 



We begin with the following four dimensional action describing gravity minimally coupled 
to a Maxwell field and charged scalar: 

C = R+-^- ^F^'^F,,, - IW - iqAi;\^ - VM) . (2.1) 

As usual we are writing F = dA, the cosmological constant is — and m,q are the 
mass and charge of the scalar field. We are interested in plane symmetric solutions, so we 
set 

ds'^ = -g{r)e-^^''^dt^ + ^ + r'^{dx^ + dy^) (2-2 

9{r) 

A = Ao{r) dt, i) = ilj{r) (2.3 

We can choose a gauge in which ip is real and work in units with L = 1. The equations of 
motion are: 

9 ^ r J 2g 



° I ^ r j ° " ^ 

X + riP'^ + ° = 2.6 

9^ 

,.U'--^),+ ':^-3r+'^^, (2,7 



r 

These equations are invariant under a scaling symmetries: 

,2, 



ar , {t,x,y) {t,x,y)/a, g a g , Aq aAo (2.^ 

6^e^, ^o^^o/^ (2.! 



Once a solution is found, this symmetry can be used to set x = at the boundary at 
infinity, so the metric takes the standard AdS form asymptotically. 
At large radius 

^o = M-f i^ = ^ + ^- (2.10) 

where A = (3 + yj9 + Am?)/2. In the boundary CFT, H is the chemical potential, p is the 
charge density, and A is the scaling dimension of the operator dual to ip. We want this 
operator to condense without being sourced, so we are only interested in solutions where 
■0 is normalizable. This typically requires setting 0(3~^) = 0. As the boundary chemical 
potential is increased beyond certain critical value, tp condenses. One may ask about the 
zero temperature limit of such a configuration. At T = 0, a condensation of il: is possible 
only if - 2q^ < -3/2. 
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2.1 = 

Here we would rephrase the results of [7] in our terms. We like to find the superconducting 
ground state of the system with a non-zero condensate. We will confine ourselves to = 
case. Being a single state without any degeneracy, a superconducting ground state does 
not have any entropy associated with it ([7], [8, 9]). 

We start by guessing a near horizon ansatz, g(r) = r'^^ip = ipQ. We use the scaling 
symmetries Eq. (2.8) to set the co-efficient in front of g to unity. Once a suitable g{r) is 
chosen, we set up a step by step perturbation in ^o- The first step is to solve for the eqn. 
of motion of Aq in this metric, 

Ao = r2+°, qi,o = (2.11) 



XI = 4^^e^v2(i+") (2.13) 



Here we have used the scaling symmetries to rescale the co-efficient of to 1. All the 
other metric component and scalar field are kept r independent at this step. 

In the next order in Aq one may solve for the r dependence of other fields (assuming 
a > —1). This procedure works as long as the various Aq dependent quantities appearing 
in the perturbative expansion are small. We get, 

ip = 'ipo-i^i{r), X = -Xi{r), g = r^-gi{r) (2.12) 

where, 



~ ~4(a + l 

9^- 4 ^ ' '^^"2(2a2 + 7a + 5) V 2 J ' ^^''^^ 

This scaling solution is valid in the regime r <^ 1. This solution may be used as a boundary 
condition to the EOM's for a numerical integration to infinity. In general the value of "0(00) 
will be non-zero. The value of a is determined by the requirement 'ip{oo) is zero. 

2.2 Small non-zero T 

From the scaling relations Eq. (2.14) one finds that IR geometry (r — )• 0) is an emergent 
AdS'^ with the same cosmological constant as that of the boundary AdS"^. At small non- 
zero T we guess that a black hole horizon will be created in the deep IR region of the 
emergent AdS^. Hence at the first step we choose the following ansatz, 

g{r) = r^(^l-^y = V'o (2.15) 

We follow the same chain of logic as in the zero temperature case and construct a 
solution as a perturbation in j4o(r). Our idea is to find out a solution such that it approaches 
the scaling solution for ^ ^ 1. As we argue below, it would then be meaningful to match 
with the numerical solution to get a full solution of the EOMs. Other quantities then 
automatically matches with the scaling solution. 
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The solution for Aq in the above metric (eqn (2.15)) is given by, 



Ao = rl+^F{-) (2.16) 
ro 



"^^"^ - r(|)r(f + f)sin(^)r^"' l^"^^ r 3 ' 3 + 3 ' 3'^ j J ('•'^^ 

Here, Aq vanishes hnearly near the black hole horizon. For r — )■ oo, F{r) ~ r^''~" + 0(r^^"'^). 
Prom that we get AQ{r) ^ r^+"(l + roO(^)) for r » tq. 

In a intermediate region r = r* ^ ro, the solutions in Eq. (2.17) could be matched 
with the zero temperature solution of ^o- Now if ro <^ 1, then we can choose the matching 
region such that ro <C r* <C 1. Importantly ^o remains small in the matching region and 
other fields could be solved perturbatively in terms of ^o (see appendix B). We get, 

Mr) = -q'i^o f'-^dr' (2.18) 



Xx{r) = -q^^lf T^df (2.19) 
Jtq 9 

51 (r) = — / r'^^dr. (2.20) 

Using the asymptotic expansion of the hypergeometric function one finds that the above 
quantities approach their zero temperature values in the matching region r ~ r* . 

Now the solution at non-zero temperature could be integrated out to infinity. As our 
solution matches with the zero temperature solution in the leading order in ^ ^ 1, an 
almost same numerical solution may be used to extend our solution to all values of r. We 
assume that the EOMs are numerically stable in the region (r,,,, oo) in a sense that a small 
perturbative change in the initial condition gives rise to a small change at infinity. 

2.3 Some results 

We do some simple calculation using our finite temperature solution. At T — )■ the horizon 
behaves like a black hole horizon situated deep inside a IR AdS"^ and various AdS"^ results 
are applicable at the leading order in ro. The temperature T of the black hole is given by, 

4^r= b'(^7exp(-xo))']^|r=ro ~3exp(-xo/2)ro, for ro < 1. (2.21) 

Hence T oc ro for small ro. Total entropy(5)^ and mass(M) of the non-superconducting 
part (i.e. for the black hole) vary as, 

5 oc rg oc (2.22) 

and, 

M oc r^ oc T^. (2.23) 



^Condensate does not have any entropy. Hence the entropy of the whole solution is same as the entropy 
of the non-supercondcuting part. 
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We may define two kind of specific lieat for our system. One is at fixed chemical 
potential(C^) and other is at fixed total charge(Cp) [24]. Here we will calculate 

dS 

C^^T—, Affixed. (2.24) 

Say a small change 5rQ in ro changes the system from (T, /^) to (T + (5T, /i + 5ii) . The 
resulting system is equivalent to a system (^if^,/^)- Here 5fi oc 6rQ and 5T oc 5rQ. Hence, 



oc oc (2.25) 
The charge density of the non-superconducting part behaves like, 

/9 0cr2^' ocrjj+° ocT3+". (2.26) 

2.4 Conductivity 

To obtain the conductivity in a background, one solves for a linearized perturbation of the 
vector potential in the same geometry. Assuming = a{r) exp{iujt) we get in our 
case. 



We define a new variable, 



In terms of this new variable. 



X 

df = —dr. (2.28) 
9 



■—Ta + V(f)a = oo^a (2.29) 



where, 



y(f ) = ^[^(2 ^ 2g2^2 exp(-x)] (2.30) 



This potential vanishes near r = 0. The superconducting nature of the system is argued 
from the existence of a supercurrent solution. If we set w = and integrate Ax from the 
horizon (with a regularity condition at the horizon), we are expected to get a non-trivial 
Ax- Existence of such a solution implies a 6 function for the real part of conductivity at 
a; = [7, 25, 26]. 

The field a has the following asymptotic behaviours near the horizon (f — )• — oo) and 
the boundary (f — )• 0): 

a{f ^ 0) al + a\f (2.31) 
a(f ^ oo) = aoe*'^^ (2.32) 
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Here, we have chosen the incoming boundary condition near the horizon. Conductivity is 
defined as follows, 

a = \ (2.33) 

It has been argued that at zero temperature, the zero frequency limit (i.e. the non- 
superconducting part) of Re((T) vanishes as powerlaw, i.e. [7, 10], 

Re((7(w)) ~ J, for 5 < 1. (2.34) 

Here 5 = V 4Vo + 1 — 1. Where Vq = limf:_5.oo r'^V{r). As the non-superconducting contri- 
bution to Re(c7) is non-zero even at small frequencies, the system does not have a energy 
gap in this channel. However non-superconducting part of Re((7(a;)) vanishes at the zero 
frequency limit. 

2.4.1 Small non-zero T 

Non-superconducting part of lim^_5.o Re((T(ci;)) is non-zero at any finite temperature. Due 
to the gaplessness of the system we expect a powerlaw decay of the above quantity with the 
temperature. The non-superconducting contribution to limi^_j.o Re(cT(cij)) has a smooth zero 
frequency limit (see appendix A) and could be calculated by setting u; = in Eq. (2.29), 
i.e. 

%a = g {A'{xf + g^e"^^^) a (2.35) 



and we have, 



(2.36) 



lim Re{a{uj)) = (2.37) 



b 



Where, and is the value of a at the horizon and boundary respectively. We break 
down the domain of r in two parts (ro,r*) and (r*,oo), s.t., tq <C <C 1. We have, 



In terms of the new co-ordinate = f. Our goal is to fix r* and take tq to a zero. In 
that case, 

2 2 

lim Re{a{uj)) = ^ ~ -rh^Ci (2.39) 

Where Ci is the limiting value of the quantity '^^'^'^ as vq — )■ 0. This value may be 
calculated from the numerics. Leading dependence of non-superconducting part of Re((j(0)) 
on ^ comes from the behaviour of the solution between (rg, r*). Taking r* in our matching 
region we can use our analytic solution in the matching region. Defining a rescaled variable 
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ri = — and the corresponding rescaled variable fi = frg, we write the equations Eq. (2.29) 



ro 

as, 



drl 



= rl{l - — ) {c^A'ix)'^ + Ig'^e'^i;^) a (2.40) 
« 2q\l{l - ^)e-^Voa (2-41) 



where we have kept the leading order terms in tq. The regular solution at the horizon has 
the following form, 



a = Im 



'0. 

Using the asymptotic expansion of the hypergeometric function we get. 



lim Re{a{uj)) ~ r^+° (2.43) 

uj—^O 



3 Non-abelian case 



The Einstein- YM action for a non-abelian gauge field with a negative cosmological constant 
is given by [5], 



C = l d'x^g [n + f^- \f^''F;)j , (3.1) 



where F^j^^ is the field strength of an SU (2) gauge field. The fully backreacted solution of 
the above equations is constructed in [14, 15]. 
The ansatz for the gauge fields is^, 

A = A{r)T^dt + B{r)T^dx. (3.2) 

To tally with the anisotropy of the gauge field ansatz in the spatial direction, we choose 
the following ansatz for our metric, 

g{r) \ J 
The Maxwell's equations of A{r),B{r) are 

Al^ yl" + A'(f + ^ + ^)-$fA = 0, 

Al^B" + B'{a^-^-^) + ^-^B = Q. (3.4) 



^Due to a repulsive term coming from the non-abelian interactions, it is expected that a isotropic ansatz 
will have a quartic instability and would possibly have more free energy than the anisotropic ones [4, 27]. 
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The diagonal Einstein equations give, 

- ' (- ^\-fl ^ ^\ 3 - ^^'2 ^^/2 gX^jl!^ 
^ \r 2c J ^yr^ rc c J 4 4r2c Ar'^gc ' 

r c \ g J g'^r'^c^ ' 

The above equations are invariant under the following scaling symmetries: 

r ^ air, {t,x,y) ^ {t,x,y)/ai, g ^ ajg, A^aiA, B ^ aiB, (3.6) 
e^^a^e^, t ^ a2t, A^A/a2- 
x^x/as, B ^ a^B. c ^ a^c. 

The second scaling symmetry may be used to set x = at infinity and the third scaling 
symmetry may be used to set c = 1 at infinity, so that the asymptotic metric is that of 
AdS^. 

The fields have the following asymptotic behavior: 

A = fi-^, B = Bl + ^, (3.7) 

where ^ is the chemical potential and p is the charge density in the boundary theory. In 
what follows we will only consider the solutions for the field B which vanishes near the 
boundary, i.e. Bq = 0. 

3.1 Zero temperature solution 

Like the abelian case the zero temperature solution is constructed by similar techniques 
[14]. We start by guessing near horizon ansatz g = r^,B = Bq, Putting this in Eq. (3.4), 
we get find out the equation of motion for A, 

r\r'A'y=^^^A = e-^, ^ = qBo/co, (3.8) 

where we have used the observation ^4 — t- at the horizon and by rescaling (3.6) we set the 
coefficients = 1- In the next order in the perturbation Aq, we get, 

B Bq - Bi{r), XG-Xi{r), 9 ^ r"^ + gi{r), c~co + ci(r). (3.9) 

All the terms with subscript 1 are sub-leading and go to zero, faster than the leading part 
where it is applicable, as r — )• 0. Here, 

XI 2r ' 4r 

where by rescaling (3.6) we may set xo = 0,co = 1. After one solves this equation by 
numerics, one again uses the rescalings of g, c, x fo make the Asymptotic geometry the 
same as that of AdSi^. For a given one numerically choose (3 in such a fashion that B 
vanishes near the boundary [14]. 
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3.2 Small non-zero T 

We follow the same strategy as in the abelian case and choose a finite temperature metric 
and B field like ^, 

5(r)=r2fl-4V B = B, (3.11) 



This satisfies the Einstein's equations at the zero'th order. Using this background fields 
we may write down the equation for A, 

+ 3, ^ = 0, /3 = ^ (3.12) 

r r*(l — rg/r'^j cq 

We would like to find a solution to the above equation which is regular at the horizon and 
approaches the zero temperature solution exp(— ^) for r ^ tq. Unfortunately there seems 
to be no analytic solutions to the above equation. However, any solution A{r) may be 
written linear combination 

A{r) = Ci exp(-^) + Ca exp(^) for r » tq. (3.13) 
r r 

We need to show that as tq — >• 0, we may choose r^, (where r* S> r^) in such a way that 
C2/C1 ^ exp(2^). Just like the abelian case we may choose a r*, s.t. ro ^ ^ 1. 
So that the nonzero-temperature solution approaches the zero temperature solution in the 
regime r ~ r.,,. This amounts to saying that Eq. (3.12) has a smooth zero temperature 
limit, at least for the solutions which are regular at the horizon. This is argued using 
matched asymptotic expansion. Near the horizon r ~ tq + 5r, (5r <C 1 and it is possible to 
linearise Eq. (3.12). The solution of the linearised equation which is regular at the horizon 
is given by, 



3/2 



2 



\ 



For small enough tq there is a region where both the above linear approximation and WKB 

3/2 

solution of Eq. (3.12) are both valid. Moreover one may choose ^ 1 in such a region. 
Using the asymptotic expansion of the Bessel function one argues that the regular solution 
matches with the correct WKB solution. Extrapolating the correct solution to r ^ tq one 
gets, 

A(r) Ci exp(— — ) for r » tq. (3.15) 

r 

This guarantees a matching region where the non-zero temperature solution approaches the 
zero temperature solution. Other fields may be solved following the similar procedure to 
that of abelian case. Using the similar argument of numerical stability we expect that our 
solution may be integrated out to infinity using a slight variation of the zero temperature 
numerics. 



^In a similar discussion of a non-zero temperature solution, the metric used was not correct in [14] and 
other formulas are also schematic. 
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3.3 Some results 



We will make some simple calculation using our scaling solution. As the black hole becomes 
a small black hole situated deep inside a IR AdS^ various AdS^ results are applicable at 
the leading order in tq. The temperature T, entropy, mass and specific heat follow the 
similar behaviour to that of abelian case. The charge density of the non-superconducting 
part behaves like, 

p cx VqA' cx exp(-— ) ~ exp(-|;). (3.16) 
3.4 Conductivity and energy gap 

In order to calculate the conductivity of this system, we need to turn on a small perturbation 
in the vector potential. We turn on the gauge field perturbations of the form: 

Al = ea{r)e~''^\'^dy (3.17) 



Here we get, 

a" + aW ^ - ^ + - + a - -^-^ - = 0. 3.18 

This can be written as a Schrodinger equation: 

- a" + V{f)a = c^Lv^a, (3.19) 

where 

,-2 132 

V{r 



g(^c'A'' + e-^^y (3.20) 



and all the derivatives in Eq. (3.19) are in terms of the new variable new variable f ( "tortoise 
coordinate") given by: 

In terms of the new co-ordinate horizon is mapped to f = — cx3. 

Here, c approaches to unity as r — >• cx), so that the spacetime is asymptotically AdS4^. 
It follows then from Eq. (3.7) that the potential V{r) vanishes near the boundary. If we 
require g ^ r"^ near the horizon at r = then the first term vanishes, while the second term 
is finite as i?(r = 0) = i?o 7^ and x is also finite at the horizon. Note that since c — )• 1 
near the boundary, the quantity w can be interpreted as the frequency of the incoming 
wave. 

Similar to abelian case, the superconducting nature of the system is argued from the 
existence of a supercurrent solution. 

The nature of the finite part of the conductivity can be inferred from the potential 
V{r) and shows a hard gap [14] at T = 0. The fact that the potential is nonzero at the 
horizon at T = makes it possible for this system to exhibit a hard gap. From Eq. (3.19), 
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the field a has the following asymptotic behaviours near the horizon (f — t- — oo) and the 
boundary (f — )• 0): 

a{r ^ 0) ~ + a\r (3.22) 
a{f ^ oo) = ooe''^^ (3.23) 



where tD = \l c^uS^ — Vq, with c/j, o/j being the near-horizon values of c and a respectively. 
Here, we have chosen the incoming boundary condition near the horizon. The conductivity 
of the system is defined in the same way as in Eq. (2.33). 
It follows from Eq. (3.19) that : 

a* a" - aa*" = 0, (3.24) 
which implies that the quantity A = a*a' — aa*' = 2ilm(aa*') is a constant. 

3.4.1 T = case 

At T=0 , equating the values of A near the horizon and the boundary we get: 

- K\' ' ^ ^ (3.25) 

0, < 

Therefore, the real part of the conductivity will vanish whenever Cj is imaginary, i.e. when 
a; < Ai = \/T^/co, which defines the gap. 

3.4.2 Small nonzero T 

However at any finite temperature T, the system does not show a hard gap as the po- 
tential V{r) actually vanishes near the black hole horizon. This is an expected behaviour 
considering thermal excitations of the condensate. One may define another gap by the 
low temperature behaviour of the conductivity. Considering the finite part (i.e. non- 
superconducting contribution) of Re(cr) at zero frequency limit at low temperature, one 
expects 

lim Re(cj(u;)) ~ exp(--f ). (3.26) 

aj-S>0 T 

Generically A2 7^ Ai and their ratio gives information about the pairing mechanism. A2 
may be thought as the mass of the charged quassiparticle carriers in the system. Ai may 
be thought as the mass of the "pairs//combination" of the quassiparticles which gives 
excitation over the pure condensate. In BCS theory such a combination of basic carriers 
is a 'Cooper pair'. In the BCS theory, Ai = 2A2. We would like to calculate A2 from our 
low temperature solution. 

The non-superconducting contribution to Re((T)(0) has a smooth zero frequency limit 
and could be calculated by setting a; = in Eq. (3.19), i.e. 

^a = ,(cV(x)2 + ,Vx^)a (3.27) 

(3.28) 



- 12 - 



and we have, 

hm Refo-ftj)) = -|. (3.29) 

Where, ah and Of, is the value of a at the horizon and boundary respectively. We break 
down the domain of r in two parts (ro,r*) and (r^,,oo), s.t., tq ^ r* ^ 1. We have, 

hm Re(a(a;)) = 4 = (3-30) 

In terms of the new co-ordinate r^, = f. Our goal is to fix r* and take tq to a zero. In 
that case, 

lim Re(.(.)) = I . -^C, (3.31) 

Where Ci is the limiting value of the quantity ""^^'2 as ro — t- 0. This value may be 
calculated from the numerics. Leading dependence of non-superconducting part of Re((T(0)) 
on ^ comes from the behaviour of the solution between (ro, r*). Taking in our matching 
region we can use our analytic solution in the matching region. Defining a rescaled variable 
ri = ^ and the corresponding rescaled variable fi = fro, we write the equations Eq. (3.27) 
as, 

-|a = -|rf (1 - ) {c^^i.f + ee--^) a (3.32) 

-^I'^^rli^-^s)--'''^- (3-33) 
'0 '1 '1 

where we have kept the leading order terms in tq. For ro ^ 1 the above equation may be 
solved using WKB approximation. One may question of validity of WKB approximation 
as the potential vanishes near the black hole horizon. This turns out not be a problem as 
we can again break down the range of ri of into two parts [1, r2] and [r2, r*). We choose our 
r2, s.t. at r ~ r2 WKB solution is valid. For a small ro, r2 lies very close to horizon ri = 1. 
For example one may choose r2 = 1 -|- ^/tq. We can also use near horizon linearisation of 
metric and solve a(ri) terms of Bessel functions for ri < r2. 



ain) = lo (^^^Efllv^ j (3.34) 

At r ~ r2 both the Bessel function and WKB method is valid and we can match these 
two following the principle of matched asymptotic expansion. Using this method one finds 
out that in the leading order only the WKB contribution matters, 

lim Re(o-(a;)) ~ exp I -2^ / '^^ I (3.35) 

'^^o \ coro Ji , ' ~ ' 
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In the above we also take — — )■ oo and — — )• 1 limit. 

ro ri 

Using the formula for temperature we get, 



3 /'°° dr 3r ( -] 

A2 = Ai— / == = Ai ;y ^ 0.669Ai (3.36) 



This may be contrasted with the weak coupling BCS value of A2 = ^Ai. We find around 
33% deviation from the BCS value. Interestingly the numerical factor is close to |. This 
gives us information about possible underlying strong coupling pairing mechanism [2]. 
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A Conductivity at u; = 

Let's consider a equation similar to Eq. (2.29), 

(P 

— ^a + V(r)a = uj'^a (A.l) 

with r E [0, 00] and the potential vanishes near the horizon (r = 00)^ The claim is 
that, 

lim Rea(w) = t^^^ (A.2) 

uj^o ^ ' |a(0)2| ^ ' 

where we find a{r) by solving Eq. (A.l) with lo = Q and using the regularity boundary 
condition at the horizon. 

From Eq. (A.l), field a has the following behaviour near boundary (r — t- 0): 

a(f ^ 0) ~ a(0) + a'(0)f (A.3) 

(A.4) 

The conductivity of the system is defined as, 

f ^ "'(0) /A r^ 

Here, we have chosen the incoming boundary condition near the horizon. 
It follows from Eq. (A.l) that : 

a*a" -aa*" = 0,, (A.6) 



*Here r is same as co-ordinate f in the main text, f is not used in the appendix due to notational 
simplicity. 
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which imphes that the quantity A = a*a' — aa* = 2ilm(aa* ) is a constant. 

Equating the values of A near the boundary and any intermediate distance rl we get: 

Re(.,.)) = M^llr™ (A.7) 

Now, for small enough uj one may break down the range of r into two parts [0, r*] , [r* , oo] . 
Where is such that V{r) S> (jj"^ for r < and uj could be treated as a small perturbation 
in this regime. Also, ^ 1 such that only leading terms contributes in V{r) for r > r^. 
In the second region we use, 

y(r)«Ciexp(-C2r) (A.8) 

We find. 

Where we have chosen the solution with correct incoming condition near r = oo. This 
solution approaches the regular solution for cj = 0, Iq ^ ^v^Ve^ — 2 ^ near r ~ r*. Also as 
3> 1, we may use the asymptotic form of the Bessel function to get, 

a(r)~e^'^'^, r~r,. (A.IO) 

Hence conductivity is given by. 



| Q(n)^ ' 
|a(0)2 



Our assertion Eq. (A. 2) is proved by noticing that as w — )• , r* — t- 00. 
B Perturbation in 

At a more technical level one may ask why we neglect the back reaction of Eq. (2.20) on the 
EOM's, i.e. what exactly we mean by Aq small. For example to neglect the back reaction 
of x'l one must have Xi^ ^ 1-) 

-r ^^^^ _ ^^■^) 

-r a (^3_i)2' i^-^i 
_ ^2+2aQ^ ^ -[^^ where Ci is a constant. (B.3) 

Where we define a = In our limit r <^ 1 although a may be large. In the last line we 

have used the the property that a^~^" (ia-i)^ bounded function in the domain [1, oo], i.e. 

there exists Ci > such that < Ci for a € [l,oo]. Using similar techniques 

one also justify the perturbative expansion for other quantities. 
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